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Electromagnons in the multiferroic state of perovskite manganites with 

symmetric-exchange striction 
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We have investigated electrically-active magnetic excitations (electromagnons) in perovskite 
manganites with the S-type (up-up-down-down) spin structure by terahertz spectroscopy. 
Eui-iYjMnOs (0.1< x <1) and Yi_ y Lu y Mn03 (0< y <1) without magnetic /-moments, which 
host collinear sinusoidal, y4-type, cycloidal, and E-type spin orders, are used to examine the system- 
atics of possible electromagnons. Three-peak structures (23, 35, 45 cm -1 ) of magnetic origin show 
up in the 25-type phase with little composition (y) dependence of frequencies, making a contrast with 
the electromagnons observed in the cycloidal-spin (x < 0.8) phases. One of these electromagnon is 
ascribed to the zone-edge magnon mode based on the calculated magnon dispersions. 

PACS numbers: 75.80.+q, 76.50.+g, 75.40.Gb 



Renewed attention to multiferroics, where the electric 
and magnetic orders can coexist, has been aroused since 
the discovery of the electric-polarization flop by apply- 
ing magnetic field in a perovskite mang anite M. The 
series of perovskite i?Mn03 (i?=Tb, Dy, Eui-^Yj) have 
a helical (cycloidal) structure of Mn 3+ spins due to the 
frustration between ferromagnetic (nearest-neighbor) J\ 
and antiferromagnetic (second nearest-neighbor) J2 ex- 
change interactions on the ab plane. The magnetically 
induced ferroelectricity has been explained by the in- 
verse Dzyaloshinski-Moriya (IDM) interaction in 
which the microscopic polarization P is expressed as 
P oc &ij x (Si x Sj) for a pair of neighboring spins Si, 
Sj with the unit vector eij connecting them. The cy- 
cloidal spins accordingly produce the macroscopic (spon- 
taneous) polarization (P s ) parallel to the spiral plane 
and perpendicular to the magnetic modulation vector q m . 
In contrast, perovskite manganites i?MnC>3 with smaller 
rare-earth ionic radius (i?=Ho, Tm, Yb, Lu, Yi-^Lu^) 
exhibit a commensurate E-type spin structure: up-up- 
down-down in the ab plane j5[, hosting also the ferro- 
electricity. The ferroelectricity in the .E-type phase orig- 
inates from the symmetric-exchange term expressed as 
Si ■ Sj, which potentially produces larger P s than IDM 
mechanism Q. In fact, a P s of 4 times as large as that 
of the a6-plane cycloidal-spin phase has been reported 
for the E-type phase @, Q. Thus, a whole series of per- 
ovskite can show versatile multiferroic states, where the 
two representative but different mechanisms -symmetric 
and antisymmetric exchange strictions- are effective. 

The dynamical interplay between magnetic and elec- 
tric excitations has also been expected for multiferroic 
materials. The first observation of electric-dipole active 
magnetic excitation, termed " electromagnon (EM)" , was 
reported for TbMnC>3 as far-infrared absorption around 



20 cm -1 [3]. Successive investigations have clarified that 
the strong infrared absorption showing up only for _E"||a 
with a two-peak structure is common to the EM in the 
cycloidal-spin phases of i?MnC>3 irrespective of the ab- 
or fee-spiral plane [Iol - [l6j . The initially predicted EM, 
which corresponds to a rotary oscillation mode of the 
spiral-spin plane (i.e. a vibration of P s ), should show 
the polarization selection rule depending on the direction 
of P s (E u \\c for P s \\a, or E u \\a for P s ||c)0, 0; how- 
ever, this is in contradiction with the experimental ob- 
servation To settle the discrepancy, the symmetric 
exchange striction (SES) mechanism has been proposed 
for the EM activity [Til H3] ; a noncollinear spin structure 
like the cycloid gives rise to the electric dipole activity 
(E^||a) of the zone-edge magnon mode via the local elec- 
tric dipole 7Cij(Si ■ Sj) built-in the i- and j-th Mn sites. 
In fact, the energy of the higher- lying mode at 60 cm -1 in 
TbMnC>3 is consistent with the zone-edge mode energy 
in the ma gno n dispersion revealed by inelastic neutron 
scattering [T9l42lj . The observed feature that the higher- 
lying EM-peak energy lowers with decreasing the radius 
of A-site ion shows a good agreement with the calculated 
result based on the Heisenberg model [22j • Recently, the 
model for lower-lyin g m ode was suggested based also on 
the SES mechanism [23|]. Since the collinear spin struc- 
ture cannot give rise to any EM in the SES mechanism, 
the suppression of the EM in the A-type phase [l2l - fl4l | 
can be reasonably explained by this mechanism. As far 
as the SES mechanism is concerned, no electric activity 
is expected for any one-magnon mode for the genuine, 
i.e. perfectly collinear, E-type phase. 

In this Letter, we have systematically investigated the 
EM in the E-type spin phase on the basis of the com- 
prehensive phase diagram established for Eui_ K Y x Mn03 
(0.1< x <1) and Yi-yLuyMnOs (0< y <1) with no /- 



2 
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FIG. 1: (color online) (a) The phase diagram of 
Eui-^YzMnOa (0.1< x <1) and Yi-yLu^MnOa (0< y <1) 
reproduced from Refs. [2a. [2r| . The Im[e/i] spectra for (b) 
A-type, (c) ab-cycloidal and (d) E-type spin phase with the 
schematic spin structures. 



electron magnetic moment on the perovskite A-site Q , by 
using terahertz time-domain spectroscopy (THz-TDS). 
In these materials, we can avoid the effect of the crystal 
field excitations of rare-earth ions, such as R —Ho, Er, 
Tm and Yb to unambiguously confirm the existence of 
the well-defined EMs in the £7- type phase. 

The high-quality polycrystallinc samples were synthe- 
sized under high pressure; the details for synthesis are de- 
scribed elsewhere || . Specularly polished samples with 3 
mm in diameter and 500 /jm in thickness were prepared. 
The terahertz wave was generated by the optical recti- 
fication of an ultrashort laser pulse with a ZnTe(llO) 
crystal and detected by a dipole antenna (24|. From 
the measured complex transmittance spectrum at 13 - 
80 cm -1 , we determined the complex optical constant 
n 2 = e^i, where the e and /i are complex dielectric and 
magnetic permittivity, respectively. In the analysis, we 
assume that \i m 1 and its negligible effect on the Fres- 
nel reflection coefficient. The validity of this assumption 
was checked in Ref . [ll[ , where the detailed experimental 
setup was also described. The possible spectral modifi- 
cation was carefully checked by comparing the spectra of 
the both single-crystalline and polycrystallinc samples of 
DyMn0 3 and Eu .8Y .2MnO 3 . 
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FIG. 2: (color online) (a) Temperature dependence of Im[e/i] 
spectra for LuMnOa. (b) Im[e/i] spectra of the _E-type phase 
at 4 K for Yi-yLu^MnOa (y=0.25, 0.5, 0.75, and 1) with 
vertical offset for clarity. Peak positions are indicated by the 
vertical dotted lines. 

Figure 1(a) shows the phase diagram of 
Eui^Yj-MnOs and Y^LtiyMnOs system [a l25l. [2t|. 
The representative spectra of the ^4-type (x=0.2), 
afc-cycloidal (a;=0.6) and E-type (y=0.5) states are 
displayed in Fig. l(b, c, d) with the respective schematic 
spin configurations. The plateau- like absorption bands 
are observed in the higher-T collincar sinusoidal phase 
for all the compounds (30 K for 2=0.2, 0.6, and 36 K for 
y=0.5) in addition to the tail of the higher-lying phonon 
modes; this broadened band has been assigned to the 
precursory EM-like excitation arising from the cycloidal 
spin fluctuation [l4l ]. The suppression of this excitation 
band is seen upon the transition to the A-type phase (4 
K) for £=0.2 (Fig. 1(b)). The sharp antiferromagnetic 
resonance (AFMR) with much smaller spectral intensity, 
which is active for the magnetic field component of light, 
is discerned at 20 cm -1 in the A-type phase In 
contrast, the development of two clear peaks (18 and 54 
enr 1 ) is identified in the cycloidal-spin phase (4 K) for 
x=0.6 (Fig. 1 (c)), similarly to the EM spectra for the 
cycloidal-spin phases of other i?Mn03 [22|]. For y=0.5 
with the E-type phase, a new EM-like feature is revealed 
in this study (Fig. 1(d)); three pronounced peaks appear 
at 24, 33 and 48 cm -1 , that are in a frequency region 
typical of magnetic (magnon) excitations. A relatively 
large spectral weight as developed below 100 K suggests 
their EM nature. 

Figure 2(a) shows the temperature dependence of 
Im[e/x] spectra for LuMnC>3 (j/=l). At 90 K (param- 
agnetic phase), only the contribution from the low en- 
ergy tail of the higher-lying phonon modes can be seen. 
The broad absorption band develops from slightly-above 
T/v=40 K with decreasing temperature, and then it shows 
rapid enhancement below T/y. This material undergoes 
the transition from the collinear sinusoidal to the £7- type 
phase at T c = 35 K with lock of q m from incommensurate 
(IC) to commensurate, accompanied by the emergence of 
ferroelectricity. In accordance with this phase transition, 
the broad absorption band turns into peak structures at 
24, 38 and 46 cm -1 with slight hardening of the frequen- 
cies while decreasing temperature. The spectra at the 
lowest temperature (4 K) for the £7- type phase are dis- 
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FIG. 3: (color online) Temperature dependence of spectral 
weight (AiVeff; see text for definition) of the electromagnon 
(or electric-dipole active magnetic excitation); (a) A-type 
(x=0.2), ab-cycloidal (z=0.6) for Eui- X YxMn0 3 , and (b) E- 
type spin phase, t/=0.25, 0.5, 0.75, and 1 for Yi-yLuyMnOs. 



played for all the measured samples (Yi_j,Lu y Mn03) in 
Fig. 2(b). All the spectra show the three-peak feature 
in common, being distinct from the feature of the two 
broad peaks in the IC cycloidal-spin phase (Fig. 1(c)). 

For more quantitative discussion, we evaluate the in- 
tegrated spectral weight per Mn site, as given by 
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where too, e and V are the bare electron mass, charge 
and the unit-cell volume, respectively. Spectral weight 
of the EM or electric-dipole active magnetic excitations 
[13] (AiVeff) is defined as the increase of N c g (between 
Wi = 13 and o>2=70 enr 1 ) from 100 K so as to eliminate 
the non-magnetic contributions. Figure 3(a) shows the 
dependence of AiVeff on the temperature normalized by 
T c =21 K for x=0.6 and by T^ FM =28 K for x=0.2 in 
Eui_ x Yj;Mn03. After the gradual increase of AN e g in 
the collinear sinusoidal phase for the both compounds, a 
rapid increase is observed for x = 0.6, while a clear sup- 
pression for x = 0.2 occurs, below the respective transi- 
tion temperatures. (A small value of AN e g (~ 0.1 x 10~ 6 ) 
for AFMR is discerned in the A-type phase for a;=0.2.) 
In contrast, AN e g in the £7- type phase below T c keeps a 
nearly constant value for all the compositions (see Fig. 
3(b)), showing again an essential difference from the EM 
in the IC cycloidal-spin phase. The electric activity of ab- 
sorption bands in the -E-type phase is again manifested 
by a fairly large AiVeff compared to that of AFMR. 

The composition (x and y in Eui.^Y^MnOa and 
Yi_ a Lu y Mn03) dependence of the frequencies of EM 
and AFMR is plotted in Fig. 4(a). In the cycloidal 
phase, the frequency of the higher-lying mode strongly 
depends on the compositions, reflecting the sensitivity of 
the zone-edge magnon energy to the A-site ionic radius 
(ta) or equivalently to the spin exchange energy [22j |. In 
fact, we can reproduce the systematic change in the fre- 
quency of the zone-edge EM mode in the light of the 
SES mechanism by assuming the exchange interactions 
as a function of ta- J a = c\a + C20TA (ol = 1,2, and 
c) in a 3-dimensional Heisenberg model (see Fig. 4 (b, 
c) for definition of J a ; c\ a and C2 a are constants) [22J. 
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FIG. 4: (color online) (a) The composition (x and y) depen- 
dence of frequencies of EM for the cycloidal and i5-type spin 
phases (or of AFMR for ^4-type spin phase). The energies of 
EM obtained by calculation are also included. The schematic 
spin configuration of the .E-type state and zone-edge magnon 
in (b) afc-plane and (c) 6c-plane. The magnetic parameters 
Ji, J 2 and J c , and unit cell length along a, b, and c axes 
are displayed, (d) The magnon dispersions calculated for the 
-E-type [9 — 0) and cycloidal phase (9 — tt/4) in the case 
of the A-site ionic radius of 1.05 A (see text). The Brillouin 
zone along is a half of that for the chemical lattice (see 
Fig. 4(b)), and hence the zone-edge mode is folded into qt~0- 
The EM is indicated by an open circle (red online), (e) Tilt- 
ing angle 9 dependence of the normalized absorption intensity 
by SES and IDM mechanisms. 



In contrast, the frequencies of all EM modes in the E- 
type phase show little dependence on ta despite the sys- 
tematic lattice-structural variation Q . In this phase, we 
propose two possible mechanisms to explain the electric 
activity of zone-edge magnon based on the symmetric- 
exchange and antisymmetric-exchange interactions. In 
the -E-type phase, the collinear spin structure parallel to 
the b axis was proposed on the basis of the neutron scat- 
tering for R=Ho, Tm, Yb and Lu 0,H^|. However, single 
ion anisotropy or DM interaction may be effective to tilt 
the spin directions off the b axis [29| , resulting in the non- 
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collincar spin state as a highly ellipsoidal cycloid. In such 
a case, the EM can show up for E u \\a via the SES mecha- 
nism as has been observed in the IC cycloidal-spin phase. 
On the other hand, the antisymmetric-exchange interac- 
tion can induce the EM in the collinear E-type state. For 
example, the zone-edge magnon mode in the collinear E- 
type ground state turns out to take the dynamical cycloid 
like spin configuration on the ab plane, which can induce 
uniform but dynamical polarization P w ||a via the IDM 
interaction, as shown in Fig. 4(b). Similarly, the zone- 
edge mode in the be plane induces P"||c, as shown in 
Fig. 4(c). Although the electric activity of these modes 
(P w ||a and P w ||c) would be distinguished by the polar- 
ization selection rule (E^Wa and -E"||c, respectively) in 
a single crystal, while it is not possible at the moment 
because of the polycrystalline nature of the samples. 

To get further insight into the nature of EM, we cal- 
culate the magnon dispersion for both the collinear E- 
type (9 = 0) and the commensurate cycloidal-spin state 
(6 = ?r/4) (see inset of Fig. 4(e) for definition of 9). As a 
typical case, we employed parameters of J\ \ — 5.6 cm -1 , 
J2 = 15 cm -1 , J c = 5.3 cm -1 and D = 2.2 cm -1 for uni- 
axial anisotropy D^(Sf) 2 for both ground states; the 
results are shown in Fig. 4(d). The zone-edge mode 
has the same frequency (~23 cm -1 ) for the both ground 
states as indicated by a circle in Fig. 4(d), resulting in 
the insensitive frequency to 9. Thus the observed lower- 
lying EM at 23 cm^ 1 can be attributed to the zone-edge 
mode for any 9. Here, the decrease of ta enhances J2, but 
hardly affects the frequency of zone-edge mode, which is 
consistent with the observed ta dependence of the EM 
frequency (Fig. 4(a)). The 9 dependence of absorption 
intensity is shown for the respective mechanisms in Fig. 
4(e). Thus, there are two possible scenarios as the ori- 



gin of EM: (i) The IDM mechanism is dominant when 
9^0. (ii) The SES mechanism, which likely gives a 
larger electric-dipole activity than the IDM mechanism, 
is dominant even for a small tilted angle 9 (pseudo-E- 
type state). The magnetic moment determined by neu- 
tron diffraction study for LuMnC>3 is 3.37 (is/Mn alog b 
axis p8| . implying the residual noncollinear component. 
Although the twin EM model [23| can be a candidate 
to induce a low-lying EM mode, too large anisotropy for 
J\ with respect to the spin direction necessary for this 
model might be unrealistic. The origin of two observed 
higher-lying modes around 35 and 46 cm -1 would be as- 
signed to the higher- lying qb=Q modes if other magnetic 
terms were taken into account. 

To conclude, we have investigated the electromagnon 
(EM) excitation in the -E-type or pseudo- -E-type (com- 
mensurate elliptical-spin) phase of the multifcrroic per- 
ovskite manganite i?Mn03 with use of non-magnetic 
A-site ions systems, Eui-^Y^MnOa (0.1< x <1) and 
Yi_ y Lu y Mn0 3 (0< y <1), by THz-TDS. The infrared 
absorption with three peaks (23, 35 and 45 cm" 1 ) is as- 
signed to the EMs in the (pscudo-) E-type phase. The 
systematic investigation on the peak positions and spec- 
tral weight revealed the distinct nature of the EM in 
the E-type phase. We proposed two possibilities for the 
electric active zone-edge magnon in the E-type phase; 
the antisymmetric-exchange (IDM) and/or symmetric- 
exchange striction (SES) mechanism. The calculation 
of the magnon dispersion demonstrates the lower-lying 
zone-edge modes, whose energy coincides with the ob- 
served lowest lying mode. 

The authors thank N. Nagaosa, T. Arima, N. Kida, M. 
Mochizuki and D. Okuyama for fruitful discussions. 
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